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Abstract: Three-dimensional second-harmonic fields, sample orienta- 
tion, and susceptibility ratios of biological samples are measured using 
polarization-resolved second-harmonic generation (SHG) microscopy. 
The three-dimensional (3D) polarization is gathered by measurement of 
a series of holograms for which excitation and analyzer polarizations 
are systematically varied, and the 3D SHG field is recovered through 
numerical back propagation. Harmonophore orientation is resolved in 3D 
from a sub- set of polarization-resolved SHG holograms. We further ex- 
pand on previous approaches for the determination of susceptibility ratios, 
adding the calculation of multiple ratio values to allow intrinsic verification. 
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Image reconstruction techniques. 
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1. Introduction 

Second harmonic generation (SHG) provides rich information regarding structural proteins and 
their organization [1-3]. This structural specificity derives from the fact that SHG signal van- 
ishes when inversion symmetry is present. This non-centrosymmetric organizational require- 
ment gives rise to an intrinsic polarization dependence in SHG [4]. This polarization depen- 
dence allows for identification of the axes of symmetry and relative strengths of SH nonlinear 
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optical tensor components of harmonophores (SH generating structures), averaged over the fo- 
cal volume of signal generation [5, 6]. 

In the context of laser-scanning second harmonic generation microscopy, this polarization 
sensitivity is exploited to obtain detailed structural information about biological tissues [7]. 
Based on models of the fibril organization, relationships between SHG intensity and incident 
linear polarization angle, orientation, and susceptibility components can be extracted through 
fitting a set of fundamental excitation angles. These analyses have been used to observe the 
orientation of collagen fibers [6], cellulose fibers [8], muscle fibers [9, 10], and cartilage [11]. 
Polarization SHG analysis can also discriminate between endogenous SHG sources from bio- 
logical tissues [12] as well as probing molecular motor dynamics [13, 14]. In combination with 
detailed models, estimates of atomic level structure, such as the helix pitch angle of the dom- 
inant harmonophore bonds can be deduced [15]. Fits to the excitation polarization dependent 
SHG intensity are time-consuming because a large number of polarization excitation angles 
must be used. With clever choice of excitation angles, a significant reduction of the amount of 
data that must be captured is achieved by using phase-cycling algebra to efficiently compute 
the fiber orientation, susceptibility terms, and harmonophore bond pitch angles [16, 17]. 

Healthy functioning of many tissues relies on proper organization, which suggests the use of 
quantitative harmonophore orientation in tissues as a pathology indicator of certain diseases. 
This approach has been pursued for identifying healthy and pathological skin dermis [18]. It 
also allows for quantitative analysis of injured tendons [19]. Polarized SHG imaging is able 
to probe collagen fiber orientation in arterial walls, which may prove useful as a diagnostic of 
coronary disease [20]. The wide range of applications of polarized-SHG (PSHG) microscopy 
for use as a quantitative diagnostic tool for diseases and other pathological conditions can be 
found in Campagnola and Dong's recent review [21]. 

PSHG is extremely powerful, and allows for observation of the three dimensional organiza- 
tion of structural tissues as the laser- scanning SHG microscopy employed restricts SHG signal 
generation to a small volume near the laser focus, providing nonlinear optical sectioning. 3D 
images are obtained by sequentially scanning a tight laser focus through a volume. The serial 
acquisition of each voxel used to construct a 3D image is a slow process, preventing application 
of SHG microscopy to rapid dynamics, restricting update rates in a 3D volume to well below 
1 Hz. This presents challenges for the direct observation of 3D temporal behavior resolved 
through PSHG, such at those present for cardiac dynamics through measurement of intrinsic 
SHG signals, which could be used to study embryonic heart dynamics [22]. To address this 
deficiency, SHG holography has been recently developed [23-26] . Holography permits capture 
of 3D object information in a 2D image, drastically improving imaging speeds. As hologra- 
phy is an interferometric technique, only co-aligned polarizations will interfere. To date, SHG 
holography methods have selected one particular pair of excitation and recording polarizations. 

In this paper, we develop a polarization-resolved SHG holographic microscope. With a set 
of SHG holograms, each with a specific excitation and recording polarization, we can recon- 
struct the orientation and susceptibility ratios of label-free, endogenous SHG generating tissue 
structures in 3D. In addition, we obtain polarization-resolved second harmonic phase images 
throughout the 3D image volume. To minimize the data that must be processed, we have adapted 
a phase-cycling approach analogous to that method pioneered by Odin et al. [16, 17] for used 
with PSHG. The resultant framework allows for computation of the total SHG intensity, and 
orientation of the local harmonophore. Our modified algorithm also admits several methods for 
computing the susceptibility ratios of the tensor elements from the SH-active tissue structures 
- providing a self-consistency check for the analysis. 
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2. Theory 



SHG holography is an interferometric technique that records the amplitude and phase of a SHG 
object field and a more powerful reference field generated by a second harmonic crystal. The 
second harmonic object field is generated through second harmonic scattering by a weakly 
focused fundamental laser beam from SH-active harmonophores in a specimen. For a given 
incident fundamental field polarization, a second harmonic field is scattered with a polariza- 
tion determined through tensor algebra that takes into account the incident fundamental field 
polarization, and the organization and orientation of the SHG harmonophore emitters in the 
specimen. Second harmonic hologram image formation is only sensitive to interference terms 
that occur through co-polarized components of the reference and object beams. Polarization- 
selection of interference presents an opportunity to fully resolve the polarization state of the 
SHG object field in three dimensions through recording and reconstruction of the amplitude 
and phase of a set of polarized-SHG holograms of a small set of 2D hologram images. 

Once the polarization state of the object field is determined, standard numerical propaga- 
tion routines provide a 3D reconstruction of the amplitude and phase of the SHG generation 
throughout the specimen volume. Further specimen information from polarized SHG data is 
possible by adopting a model of the local second-order nonlinear susceptibility tensor. As 
SHG requires media that lack inversion symmetry on both microscopic and macroscopic scales, 
SHG signals are restricted to tissues that display this organization. Prior work has studied the 
organization of many second harmonic generation active tissues. A number of groups have 
shown that local susceptibility tensor can be described in terms of hexagonal or cylindri- 
cal symmetry, as demonstrated for muscle [9, 15], collagen [27,28], and starch [29]. Using 
the contracted notation d[ s = 1/2£qX^ , the possible non-zero tensor elements are given as 

d 33 = l/2e 0 x£l d 31 = l/2e 0 ^z ( S = l/2£oXw, d ^ = V^iS = l/2£oxH = l/2e^Xy% 

du = l/2e 0 Xzyl = l/2£ 0 ^xzy = -l/2e 0 ^yxz = -l/2e 0 Xyll. Assuming Kleinman symmetry 
provides further simplification, with du = 0 and d\$ = d 3 \. The ratio of the two remaining 
tensor components, p = gives information about the ordering and supramolecular organi- 



While we are making use of femtosecond laser pulses of a few hundred femtoseconds in du- 
ration, we may theoretically describe the SHG hologram formation by considering a constant- 
wave (cw) fundamental field for illumination, E- m . This incident field is assumed to be linearly 
polarized along the x-axis of the transverse coordinate frame of the specimen. In a coordinate 
frame (x,y,z) in the region of the specimen, where z is the direction of propagation of the fields 
and assuming that the field is not localized strongly enough to give rise to a strong longitudinal 
field, we can write the SHG polarization density generated by an incident field polarized lin- 
early along an angle a with respect to the x-axis by a fiber oriented at an angle 6 again with 
respect to the x-axis as 



where c(r) is the spatial distribution of the harmonophore concentration normalized, to the 
peak value in the specimen, and r = xx + yy + zz is the position vector inside of the specimen. 
Here, we see that the SHG source term generating the SH scattered object field is proportional 
to the local normalized harmonophore concentration and local specimen orientation. These 
polarization densities lead to a second-harmonic field scattered by the incident field. We record 
the generated field as a coherent superposition of second-harmonic scattered from within the 
volume of the specimen illuminated by the fundamental field. We can record a hologram of this 



zation. 




E? n sin [2a + 20 (r)]rf 3 ic(r) 

El {sin [a + 0(r)] 2 d 3l + cos [a + 0(r)} 2 d 33 } c(r) 



(1) 



(2) 
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Fig. 1. Reconstructed field magnitude of 16 holograms of a 20 /im thick slice of corn seed, 
with varying excitation and analyzer polarizations. Rows vary excitation angle from 0 to 
135° in 45° steps, from top to bottom; columns vary analyzer angle from 0 to 135° in 45° 
steps, from left to right. The colormap indicates total SHG intensity, with low intensity in 
black and high intensity in white. 



field with the inclusion of a reference field [24] . An analyzer set at an angle \\f with respect to 
the x-axis selects the polarization of the field incident on the camera as a superposition of the x 
and y SHG fields, as 

E SHG = E SUG cos ( VO + E SHG sin ( VO ■ ( 3 ) 

We can also control the polarization of the illumination field, allowing us to set its polarization 
at some angle a, again with respect to the x-axis. We record a set of holograms with varying 
illumination and analyzer angles, a = ^ and iff = where n and m are both integer values 
between 0 and 3. Using standard holographic processing [24, 30-33] to back-propagate the 
holograms we can recover the field at any z distance. We can write each field at this distance as 
E m ,n = jc{*)F mi n exp (i<f> m ,n) exp (i(j> m ^n)- With the normalized relative field magnitude given by 



: cos (_) sin + 20J + - sin (-) [(1 + p) - (1 -p) cos ( — + 



(4) 



1 

4=2 £ l^| 2 = f/ + Vcos(^ r )+Wcos(^7T + 4e) (5) 



V 2 

we can then write the total SHG intensity in the back-propagated plane as 

=l 

with[/= \ (3p 2 + 2p+7),V = \ (p 2 -l),andW= \ (p 2 -2p-3) defined as in [17]. Solv- 
ing for the parameters in terms of measured intensities, we arrive at equations for the three 
parameters in terms of the above-measured intensities, 

tf=^(/o+/i+/2+/ 3 ), V = ±^a/a/ 3 2 1 +A/ 0 2 2 , W = I(A/oi+A/ 23 )|^| (6) 
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where AI a b —I a —l\). We can recover the sample-orientation angle 0 for each point in the recon- 
struction from these measurements as tan (26) = A/31 /A/02 • As this angle estimate is derived 
from a ratio of intensities, it is insensitive to the local concentration of the harmonophore or the 
collection efficiency of the SH light. Making use of ratios eliminates local concentration and 
SHG intensity factors, allowing for extraction of local susceptibility parameters such as p and 
6 without requiring knowledge of harmonophore concentration distributions in the specimen. 
Similarly, solutions for the value of the susceptibility ratio p can be formed from ratios of the 
any pair of the intensity terms U,V, and W. While solutions of p should be independent of the 
calculation method, noise in the data present variations in the values extracted by each algo- 
rithm. Cross-checking multiple solutions for consistency of the extracted parameters validates 
the solutions of p (r) . 

3. Experiment 

The holography setup is a modification of our previously described system [24]. A Yb:KGW 
oscillator producing 340 fs pulses centered at 1027 nm is directed into a modified Mach-Zender 
interferometer. The reference field is generated in a 100 /im thick KH2PO4 (KDP) crystal. 
The reference is separated from the fundamental using a dichroic filter, and directed to the 
camera through a spatial filter and adjustable time-delay arm. The remaining fundamental light 
passes through a half-wave plate for excitation polarization control and is focused onto the 
sample using a Meiji 4x 0.1 NA objective. The SHG signal is collected in the forward- scattered 
direction using a Zeiss Epiplan 50x 0.5 NA objective. The object and reference SHG beams are 
combined using a non-polarizing beam cube and pass through a wire-grid analyzing polarizer 
before being collected by an EM-CCD camera. This polarizer is rotated to measure different 
polarization projections on the camera. 

We used as samples prepared slides of 20 fim thick sections of corn seed and canine tongue. 
The polarization of the horizontal linearly-polarized beam incident on the sample is rotated 
through 0°, 45°, 90°, and 135°. The resulting field is then joined with the reference, and the 
reference polarizer rotated through 0°, 45°, 90°, and 135°. This results in a collection of 16 
holograms, each of which is back-propagated numerically to yield the field in the image plane. 
An example of the recovered field intensity is shown in Fig. 1, taken with a 20 fim thick slice 
of corn seed as the sample. 

From these 16 fields, we calculate the four second-harmonic intensities, I n , by summing 
the intensity of each of the recovered field for all excitation polarizations for each analyzer 
polarization. From the four analyzer intensities, we can determine the value of the parameter U . 
We can then calculate the other parameters, V and W, from differences between the recovered 
intensities. 

4. Results and discussion 

From the set of 3D reconstructed SHG fields, we can obtain 3D image distributions of the total 
SHG scattered intensity (not modulated by the local projection of the incident polarization onto 
the symmetry axis of the organized tissue), the local SHG tensor symmetry axis orientation, 0, 
the SHG field phase, and the local value of p . We will first concentrate on the first three quan- 
tities, and show 3D image data for SHG intensity, orientation, and phase at each z-plane in the 
sample. The total SHG intensity is computed from the data by /sHG,tot = \ X^=o^ = ^> w ^ m 
the local orientation computed through 6 = | tan -1 (A/31 /A/02). The local orientation is shown 
in Fig. 2(a-b), with lines superimposed on the total intensity above a SHG intensity reaching 
14.5% of the peak SHG intensity, i.e., for /sHG,tot > 0.145 x max(/sHG,tot)- The phase difference 
between the SHG fields for horizontal excitation polarization and 45° analyzer polarization and 
vertical excitation and analyzer polarization, taken at the reference back-propagation plane, is 



#170543 - $15.00 USD Received 13 Jun 2012; revised 30 Jul 2012; accepted 31 Jul 2012; published 3 Aug 2012 
(C) 2012 OSA 1 September 2012 / Vol. 3, No. 9 / BIOMEDICAL OPTICS EXPRESS 2009 



0 0.5 1 -2n -71 0 7i 2n 

Fig. 2. Frame from a movie showing reconstructed intensities of canine tongue (a, Media 
1) and corn seed (b, Media 2), with polarization orientation indicated with white lines, 
showing reconstructions over a distance of 22 jirni through the sample. Phase differences 
are shown for tongue (c) and corn seed (d), between projections (n,m) = (1,2) and (0,0). 



shown in Fig. 2(c-d) for the threshold region as described above. The back propagation distance 
for these holograms is 58 /im. We display the SHG phase difference for only one pair of the six- 
teen polarization/excitation components, with the difference taken to remove the propagation 
phase, as we have found no significant variation in the phase behavior of the 16 polarization- 
isolated SHG holography terms. This suggests a lack of polarization dependence on SHG in- 
tensities, and that no birefringent phase matching effects are apparent in the specimens that we 
examined. 

Determining the value of the susceptibility tensor ratio components is more challenging than 
the extraction of local orientation. We found solutions for p from each of the local intensity and 
concentration independent ratios of ryu = V/U, rwu = W/U, and ryw = V/W. As there is a 
sign ambiguity in V, ryu provides four solutions for p, while ryw and ryu each provide two 
solutions, leading to a total of eight solutions for p in a single data set. 



-rvu±2 A /4 + 4r V u-5r 2 



3r vu + 4 



vu 

(7) 



-r V u ± 2\ 4 - 4r V u - 5r 2 

(8) 



3rvu-4 



-(1 + rwu) ±2^/1-5^ 
^ 3 r W u - 1 

P = (10) 
ryw ±4 

The value of p is calculated for each pixel in the reconstructed SH field using the three different 
ratios (and the two possible roots, for the cases of ryu and rwu)- From this array of values, we 
calculate the mean and standard deviation of all non-zero values contained in these retrieved 
p images. The multiple p solutions to provide a self-consistency check on the values of p 
obtained. In the data that we have analyzed to date, we found that the V-W ratio, ryw, provides 
the most consistent results, perhaps because it does not require the choice of a root in the 
solution and all values remain real. The only ambiguity is the sign of V. For the tongue sample, 
as seen in Fig. 3, we determine four values for p. From the V-W ratio, with V > 0, we find 
p = 3.4, a = 1.3. From the V-W ratio, with V < 0, we determine p = 2.9, o = 1.4. We can 
use the other ratios to verify this result. From the W-U, positive root, we get p = 3.4, a = 1.3, 
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Fig. 3. Calculated p values per pixel for two example ratios: (a) V/W, V > 0 and (b) V/W, 
V < 0. The histogram of ratio V/W, V < 0 is shown in (c). All values are displayed over 
the same range. 



while the negative root of the same ratio lacks sufficient non-zero values to yield a p value. 
From the V-U ratio, V > 0 with a positive U root, we retrieve the value p = 2.2,cr = 0.35, 
while the V < 0 and both negative root solutions have so few non-zero datapoints as to not 
provide a solution. We see good agreement between the two values determined from the V-W 
ratio, which is further verified by agreement with the value determined from the W-U ratio. 
The divergent value, determined from the V-U ratio, can then be given less weight in light of 
the good agreement of the other values. 

5. Conclusion 

We have demonstrated a polarization-based SHG holographic microscope, which we used to 
gather sample orientation and susceptibility ratios in three dimensions from unlabeled biologi- 
cal samples from a set of excitation and analyzer polarization-resolved SHG hologram images. 
From these data, we have numerically reconstructed three dimensional images of the total (not 
modulated by the incident field polarization) second-harmonic intensity scattered by the spec- 
imens. The holographic data provide access to the SH scattered field phase throughout the 
specimen volume. We found no polarization dependence on the phase images, which strongly 
indicates that there are no polarization-dependent phase variations, ruling out any birefringent 
phase matching effects on the SHG hologram image formation. The 3D distribution of sample 
orientation and the susceptibility ratio, p, of the local variations in %^ tensor are also extracted 
from analysis of the PSHG hologram images. The methods presented have inherent consistency 
checks, and the agreement between the different methods has been examined as a validation 
metric. While a full analysis is presented here, it is clear that a sub-set of polarization-resolved 
SHG holograms can be recorded in order to extract 3D images of the total SHG intensity, orien- 
tation, optical phase, and relative susceptibility values. Given the ability to capture 3D images 
in a rapid time (each image could be acquired in less than a millisecond), this new PSHG holog- 
raphy opens the pathway for high-speed dynamical 3D PSHG imaging for studying biological 
processes, such as muscle mechanics with millisecond temporal resolution. 
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